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New kind of differential equations, called local fractional differential equations, has been pro- 
posed for the first time. They involve local fractional derivatives introduced recently. Such equations 
appear to be suitable to deal with phenomena taking place in fractal space and time. A local frac- 
tional analog of Fokker-Planck equation has been derived starting from the Chapman-Kolmogorov 
condition. Such an equation is solved, with a specific choice of the transition probability, and shown 
to give rise to subdiffusive behavior. 
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Derivatives and integrals of fractional order have found 
many applications in recent studies of scaling phenom- 
ena |[]-|4j . The main aim of the most of these papers is to 
formulate fractional integro-differential equations to de- 
scribe some scaling process. Modifications of equations 
governing physical processes such as diffusion equation, 
wave equation and Fokker-Planck equation have been 
suggested [^10 which incorporate fractional derivatives 
with respect to time. Recently Zaslavasky |11 argued 
that the chaotic Hamiltonian dynamics of particles can 
be described by using fractional generalization of the 
Fokker-Planck-Kolmogorov (FPK) equation. However 
fractional derivatives are nonlocal and hence such equa- 
tions are not suitable to study local scaling behavior. In 
the present work we rigorously derive fractional analogs 
of equations like the FPK equation involving one space 
variable. Our approach differs from the above mentioned 
ones since we use local fractional Taylor expansion, which 
was established only recently |l2[ . As is argued below, 
such equations can provide appropriate schemes for de- 
scribing evolutions (e.g. sub or superdiffusive) normally 
not obtained from the usual FPK equation. 

It was realized recently [|l2| that there is a direct 
quantitative connection between fractional differentiabil- 
ity properties of continuous but nowhere differentiable 
functions and the dimensions of their graphs. In order 
to show this, a new notion of local fractional derivative 
(LFD) was introduced. The LFD of order q of a function 
f(y) was defined by 

Wf(y) = lun d yM- f M ] 0< q <l (1) 
x^y [d{x — y)\ q 

where the derivative on the RHS is the Riemann-Liouville 
fractional derivative viz., for < q < 1 



d«f(x) 



[d(x-a)]i T(l-q)dxJ a {x-y)i 
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As is obvious from equation (g), the opera- 
tor d q /[d(x — a)] q is nonlocal and further the 
d q f(x)/[d(x-a)] q ^ for f(x) = Const. The moti- 
vation for the definition of E) q was to correct for both of 



these features. It was shown in [fl2f , in particular, that 
the LFD of Weierstrass nowhere differentiable function 
exists upto (critical) order I— 7, where I+7 (0 < 7 < 1) is 
the box dimension of the graph of the function. Further 
the use of LFD to study pointwise behavior of multi- 
fractal functions was also demonstrated. The definition 
was generalized Jl5| for a function for which the first N 
derivatives exist by replacing \f(x) — f{y)\ on RHS of 
equation ([!]) by 



F N (x,y) = f(x) - 



N 

n=0 v 
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with q in the interval (N,N + 1]. Sometimes it is essen- 
tial to distinguish between limits taken from above and 
below. In that case we define 



E>if(y)= Urn 



d q F N (x,y) 
y±[d±(x-y)]i' 



(4) 



We will assume E) q = 1D\ unless mentioned otherwise. 
We note that when q = n, an integer, K) q f(y) is simply 

'i". ran '///"• 

The importance of the above definition also lies in the 
fact that the LFDs appear naturally in the fractional 
Taylor expansion as the coefficient of the power with frac- 
tional exponent. Thus, for A = x — y 



N 



f {n) (y) 

T(n + 1)' 



]D q ±f(y) 

T(q + 1) 



(±Ay + R q (y,A) (5) 



where R q (y,A) is the remainder p2[| . 

The basic idea of the present paper is to utilize 
such fractional Taylor expansions in the Chapman- 
Kolmogorov condition and obtain analogs of FPK equa- 
tion. We begin by recalling the usual procedure and dif- 
ficulties of obtaining FPK equation. Let W(x,t) denote 
the probability density for a random variable X taking 
value x at time t, then 



W(x,t 



P(x,t + t\x' \t)W{x' ,t)dx' 



(0) 
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where P(xi, ti\x2, ti) denotes the transition probability 
from x\ at time t\ to X2 at time ti and r > 0. The 
usual FPK equation is obtained from equation (||) 
by expanding the integrand in a Taylor series. 

There are number of limitations of this approach aris- 
ing naturally from the assumptions going into its deriva- 
tion. For instance, as noted in p7| , probability distri- 
butions whose second moment does not exist are not 
described by FPK equation even though such distribu- 
tions may satisfy original Chapman-Kolmogorov equa- 
tion. Also, as emphasized in ]l8[ |, the differentiability as- 
sumption may also break down in various situations. For 
instance, the transitional probability density may not be 
differentiable at x = x' in which case the derivation of 



FP equation itself will break down. Another situation is 
when we have a fractal function as the initial probabil- 
ity density. In such a case even the usual Fokker-Planck 
operator can not be operated on the initial density. 

It is thus of interest to broaden the class of differential 
equations one can derive starting from the Chapman- 
Kolmogorov equation, and study various processes de- 
scribed by them. In this paper we pursue the possibility 
of removing the assumption of differentiability of prob- 
ability densities. We follow the usual procedure to de- 
rive the Fokker-Planck equation from equation (^|) except 
that we now expand the integrand using fractional Tay- 
lor expansion (^|) instead of ordinary Taylor expansion. 
Thus, if A = x — x' , 



W(x,t + - 



= W(x, t) + Y / W^-T,(7^—X I dx'A»P(x + A,t + t\x, t)W(x, t) 
^r(n + l) d(-x) J 

dy(y-xfP(y,t + T\x,t)W(x,t)] 

dy(x — yf P(y, t + t\x, t)W(x, t)] + Remainder 



1 
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where E) x is a partial LFD w.r.t x. Now if < a < 1 

W(x, t + r)- W(x, t) = y K^- + Remainder. 

T(a + 1) 

where E) t is partial LFD w.r.t. t. In general a and (3 
may depend on x and t. But we assume that a and /3 
are constants. Therefore we get 

r (a +i) ^h d^y [ r (" +1 ) ] 



a r Mt(X,t,T) 



x+[ r(/3 + i) 



W(x,t)] (8) 



where 



M+(x,t,r)= dy(y-x) a P(y,t + T\x 1 t)a>0, (9) 



M-{x,t,r)= dy(x~y) a P(y,t + T\x 1 t)a>0 (10) 

J — OO 

and 

M a (x,t,r) = M+{x,t,r)+M-(x,t,T) (11) 

are transitional moments. The limit r — > gives us a 
equation 



E>?W(x,t) = C(x,t)W(x,t) 
where the operator C is given by 



(12) 



N 



C ^ = E (g^r) n A n a (x 1 t) + lD^AU^) 



+E>P + Ai + (x,t) 



(13) 



where 



A?(x, t) = Um P ' \ (14) 



and 



Ai{x,t)=A p a+ {x,t)+Ai_{x,t). 



(15) 



Here corresponding A Q 's are assumed to exist. We would 
like to point out that the equation (|l^) is analogous to 
truncated Kramers-Moyal expansion. Two rather impor- 
tant special cases are, < (3 < 1 and 1 < (3 < 2. In the 
former case we get the operator 

C(x,t) = lD x _A^(x,t) + Z£X+0M), 

and in the latter case we get 



9 A i 



Cu.l) -^A^x,t)+n^_Ai_{x,t)+n^ + Ai + {x,t) 
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This operator can be identified as generalizations of the 
Fokker-Planck operator in one space variable. It is clear 
that when a = 1 and (3 = 2 we get back the usual Fokker- 
Planck operator. 

It may be pointed out that the local fractional differ- 
ential equations (LFDE) that we are proposing here are 
new kind of differential equations. To our knowledge this 
is the first direct occurrence of such equations. We note 
that they are different from the conventional fractional 
differential equations which have been studied to some 
extent in the literature 1 13 14J] and which have found sev- 
eral applications ranging from solutions of Bessel equa- 
tion, diffusion on curved surfaces to wave equation etc. 
In fact the equations appearing in are an con- 

ventional fractional differential equations. On the other 
hand the present LFDE involve operators -ZD 9 , which 
found successful applications Jl2| in studying differentia- 
bility properties of nowhere differentiable functions and 
relating them to dimensions. They are appropriate to 
address scaling phenomena. It is for this reason that one 
would expect the equations governing the fractal pro- 
cesses to be LFDE. At this stage it is worth reflecting 
for a moment on the behavior of meaningful solutions of 
simple LFDE. We begin by considering the equation 



E>lf(x)=g{x). 



(16) 



The questions of the general conditions guaranteeing so- 
lutions of such an equation is an involved one. We note 
that the equation E)%f{x) — Const, does not have a fi- 
nite solution when < q < 1. Interestingly, the solutions 
to (|l6|) can exist, when g(x) has a fractal support. For 
instance, when g (x) = xc(x), the membership function 
of a cantor set C (i.e. g(x) = 1 if x is in C and g(x) = 
otherwise), the solution with initial condition f(0) = 
exists if q = a = dim#C. Explicitly, generalizing the 
Riemann integration procedure, 



f(z)= Pc(x \= Hm y %i 

v ' 2=0 V ' 



■Fb 



(17) 



where subdivision points of the interval [x = 

0,xn = x] and Ffc is a flag function which takes value 1 if 
the interval [xt, Xi+i] contains a point of the set C and 
otherwise. Note that Pc(x) is a Lebesgue- Cantor (stair- 
case) function and satisfies the bounds ax a < Pc(x) < 
bx a where a and b are suitable positive constants. In 
general, the algorithm of the equation (^j) will work only 
for the sets C for which dim^C = dim^fC (in fact in this 
case only N a terms in the summation are nonzero). More 
details about solutions of such equations and algorithms 
will be discussed elsewhere [jl9| . 

Returning back to equation (|l4|) it is clear that the 
small time behavior of different transitional moments de- 
cide the order of the derivative with respect to time ( in 
order to demonstrate this point we consider the example 



of a Levy process below). On the other hand, small dis- 
tance behavior of transitional probability or the differen- 
tiability property of the initial probability density would 
dictate the order of space derivative. Depending on the 
actual values of a and (3 as well as their interrelation the 
above local FFPK equation will describe different pro- 
cesses. 

Equations which give rise to evolution-semigroup are 
of interest in physics. The equation ( |l2| ) corresponds to 
a semigroup if a = 1. One can then write down a formal 
solution of the above equation in this case as follows. In 
the time independent case we have 



W(x,t) = e £ ^*l/F(x,0) 



(18) 



and when the operator depends on time we have 

W(x, t) =? X C{x - t,)dt 'w{x, 0) (19) 

where £ is an operator in equation (|l^) and T is the time 
ordering operator. 

For the symmetric stable Levy process of index /i, the 
moments scale as M 7 (Ai) = A 7 ^M 7 (t) and we get 

1DJ^W( X , t) = 1DI_ (x, t)W(x, t)] 

+JD2 + [AJ f/l ^{x,t)W{x,t)]. (20) 

Since the process is symmetric the first derivative does 
not appear. The order of the time derivative depends 
on that of space derivative but it is always less than one. 
Now there is only one free parameter 7 which is restricted 
to the range < 7 < /i. In this case the value of 7 will 
be decided by the differentiability class of the initial dis- 
tribution function. (The details and intricacies will be 
addressed in p9[ .) When \i = 2 and 7 = 2 we get back 
the usual Fokker-Planck equation describing a Gaussian 
process. Equation (|2^) forms one example where usual 
derivation of FPK equation breaks down and we get non- 
trivial values for the orders of the derivatives. 

As our next example we consider the transition prob- 
ability 



P(x,t + r\x',t) 



1 



(21) 

0rAP c (i,r) 

= 8{x-x') if AF c (t,r) = 0, (22) 



where AP c («,t) = P c {t + r) - Pc(t))- This transition 
probability describes a nonstationary process which cor- 
responds to transitions occouring only at times which lie 
on a fractal set. Such a transition probability can be used 
to model phenomenon where transition is very rare, for 
instance, diffusion in the presence of traps. The second 
moment is given, from equation (|ll]), by 



M 2 (V 



AP C (V 



-Xc(t) 



!E> a P c (t) 
2 r(a + l) 



(23) 
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This gives us the following local fractional Fokker-Planck 
equation (in this case an analog of a diffusion equation). 



JD?W(x,t) 



—Xc(t)^W(x,t) (24) 



We note that even though the variable t is taking all real 
positive values the actual evolution takes place only for 
values of t in the fractal set C. The solution of equa- 
tion (p4|) can easily be obtained as 



W{x,t) = P t - t0 W(x,t a ) 

where 

N - 1 i n2 



(25) 



(26) 



be of value in the studies of anomalous diffusion, chaotic 
Hamiltonian systems, disordered phenomenon, etc. In 
our derivation we assumed that the orders a and f3 of 
derivatives involved are constants. This would require a 
modification the description of multiscaling multifractal 
processes. We further note that directional LFDs are de- 
fined in j20). Using them it may be possible to obtain 
local FPK equation involving several variables. 
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The above product converges because except for number 
of terms of order N a all other terms take value 1. It is 
clear that for to < t' < t 



W(x,t) = P t - V Pt'-t W{x,t ) 



(27) 



and P t gives rise to a semigroup evolution. Using equa- 
tion ([It]) it can be easily seen that 



W(x,t) = e^"^W(x,t = 0). 



(28) 



Now choosing W(x,0) = 5(x) and using the Fourier rep- 
resentation of delta function then we get the solution 



W(x,t) = 



1 



-e p c(*) 



(29) 



Its consistency can easily be checked by substituting this 
in Chapman-Kolmogorov equation. This solution satis- 
fies the bounds 



1 -* 2 1 
e w < W(x, t) < 



(30) 



for some < a < b. This is a model solution of a sub- 
diffusive behavior. It is clear that when a — 1 we get 
back the typical solution of the ordinary diffusion equa- 
tion which is (7rf) -1 / 2 exp(— x 2 /t). 

To conclude, we have derived the generalization of 
FP equation which involves the local fractional deriva- 
tives. Our equations are fundamentally different from 
any of the equations proposed previously since they in- 
volve LFDs. They are local and more natural generaliza- 
tion of ordinary differential equations. LFDEs deserve a 
separate study in their own right jl9| . We would like to 
point out that LFDEs naturally give rise to dynamical 
systems of new kind (neither discrete nor continuous) in 
which time evolution takes place for values of time be- 
longing to a Cantor- like set. We further remark that since 
the present analogue of FPK equation is derived from 
first principles, we feel that our equation will have general 
applicability in the field of physics. We expect them to 



* Electronic address: kirkol@physics.unipune.ernet.in 

* Electronics address: adg@physics.unipune.ernet.in 

[1] T. F. Nonnenmacher, J. Phys. A: Math. Gen. 23, L697 
(1990). 

[2] W. G. Glockle and T. F. Nonnenmacher, J. Stat. Phys. 

71, 741 (1993). 
[3] M. F. Schlesinger, J. Stat. Phys. 36, 639 (1984). 
[4] B. B. Mandelbrot and J. W. Van Ness, SI AM Rev. 10, 

422 (1968). 

[5] M. Giona and H. E. Roman, J. Phys. A: Math Gen. 25, 
2093 (1992) 

[6] H. E. Roman and M. Giona, J. Phys. A: Math. Gen. 25, 

2107 (1992) 
[7] W. Wyss, J. Math. Phys. 27, 2782 (1986). 
[8] W. R. Schneider and W. Wyss, J. Math. Phys. 30, 134 

(1989). 

[9] G. Jumarie, J. Math. Phys. 33, 3536 (1992). 

[10] H. C. Fogedby, Phys. Rev. Lett.73, 2517 (1994). 

[11] G. M. Zaslavsky, Physica D 76, 110 (1994). 

[12] K. M. Kolwankar an d A. D. Gangal, Chaos 6 505 (1996). 
( |chao-dyn/9609016| ) 

[13] K. B. Oldham and J. Spanier, The Fractional Calculus 
(Academic Press, New York, 1974). 

[14] K. S. Miller and B. Ross, An Introduction to the Frac- 
tional Calculus and Fractional Differential Equations 
(John Wiley, New York, 1993). 

[15] K. M. Kolwan kar and A. D. Gang al, Pramana - J. Phys. 
48, 49 (1997) . ( |chao-dyn/971 101$ 

[16] H. Risken, The Fokker-Planck Equation (Springer- 
Verlag, Berlin, 1984). 

[17] E. W. Montroll and M. F. Schlesinger, On the wonder- 
ful world of random walks in Nonequlibrium Phenomenon 
II: From Stochastic to Hydrodynamics edited by J. L. 
Lebowitz and E. W. Montroll (North-Holland, Amster- 
dam, 1984) 

[18] W. Feller, An Introduction to Probability Theory and its 

Applications (Wiley, New York, 1968) Vol 2. 
[19] K. M. Kolwankar and A. D. Gangal, (in preparation) 
[20] K. M. Kolwankar and A. D. Gangal, in the proceedings 
of the Conference 'Fra ctals in Engineering', Archanon, 
1997. (|physics/9801010D 



4 



